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Beschreibung 

Verfahren und Anordnung sowie Computerprogramm mit Programm- 
code-Mitteln und Computerprogramm-Produkt zur Ermittlung ei- 
ner ausgewahlten Position einer mobilen Kommunikationsein- 
richtung in einem Kommunikationsnetz 

Die Erfindung betrifft eine Ermittlung einer Position einer 
mobilen Kommunikationseinrichtung in einem Kommunikationsnetz 
(Lokalisierung) . 

Mit zunehmender Verbreitung mobiler Kommunikation nimmt auch 
eine Nachfrage nach zusatzlichen Dienstleistungen, sogenann- 
ten „services*, bei Mobilfunksystemen zu. 

Unter ^Location Based Services' werden dabei zusatzliche 
Dienstleistungen von Mobilfunkbetreibern verstanden, welche 
Nutzern von den Mobilfunksystemen ortsgebunden, d.h. abhangig 
von einer Position oder einem Standort des jeweiligen Nut- 
zers, angeboten bzw. bereitgestellt werden, beispielsweise 
orts- oder entf ernungsabhangige Nutzungstarif e oder Orientie- 
rungshilfen fur Rettungseinsatze oder Suchdienste. 

Grundlegend fur einen ^Location Based Service* ist demzufolge 
eine Lokalisierung oder Positionsermittlung des jeweiligen 
Nutzers bzw. dessen mobilen Kommunikationsgerats . 

Fur eine solche Lokalisierung von mobilen Kommunikationsgera- 
ten in Kommunikationsnetze sind verschiedene Techniken be- 
kannt, beispielsweise eine Positionsermittlung auf Basis ei- 
ner Lauf zeitbestimmung bzw. Lauf zeitmessung von Kommunikati- 
onssignalen von einem mobilen Kommunikationsgerat zu einer 
Basisstation e'ines Kommunikationsnetzes ( [1] , [2]) oder eine 
Lokalisierung mittels satellitengestutzter Systeme, wie ein 
GPS. 
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Das aus [2] bekannte, Laufzeit basierte Positionsermittlungs- 
verfahren wird fUr ein Mobiltelefon, im allgemeinen einer Mo- 
bilstation, in einem GSM-Kommunikationsnetz (= Global System 
for Mobile Communications) ([7], [8], [9]) gemafi einer TDMA- 
Mobilfunk-Technik (Time Division Multiple Access-Mobilfunk) 
durchgefiihrt . 

Einer einzelnen Mobilstation, die sich bei einer ortsfesten 
Basisstation (gesprachsfiihrenden Basisstation) eingebucht 
hat, wird bei dieser eine derzeit freie Zeitlage in einem 
TDMA-Rahmen zugeteilt. ^ 

Zu dieser Zeitlage gehen die ftir die betreffende Mobilstation 
bestimmten Kommunikationssignale in Signalpaketen, sogenann- 
ten bursts, mit einer Lange von 15/2 6ms von der Basisstation 
ab, beziehungsweise mtissen die von der Mobilstation gesende- 
ten Kommunikationssignale bzw. bursts an der Basisstation 
eintref fen. 

Die von der Basisstation emittierten Kommunikationssignale 
finden ihren Weg zur Mobilstation aufgrund von Streuungen tt- 
ber unterschiedliche Wege (Mehrf achausbreitungen) , wobei sie 
frequenzabhangig gedampft werden. ^ 

Eine Empf angsf eldstarke der von der Mobilstation empfangenen 
Kommunikationssignale ist somit nicht nur von einer Entfer- 
nung der Mobilstation von der Basisstation, sondern auch von 
der Frequenz und topographischen Gegebenheiten zwischen Mo- 
bilstation und Basisstation abhangig. Deshalb werden die ein- 
zelnen Datenpakete auf verschiedenen Tragerf requenzen ver- 
schickt, wodurch selektive St5rungen einer Frequenz auf meh- 
rere Teilnehmer verteilt werden kOnnen. 
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Daftir ist allerdings eine prazise Synchronisierung zwischen 
Mobilstation und Basisstation notwendig. Diese Synchronisie- 
rung wird durch eine Mobilitat eines Nutzers zusatzlich er- 
schwert, weil die Mobilstation sich nun in sich andernden 
Entfernungen zur Basisstation befindet und ihre Kommunikati- 
onssignale unterschiedliche Laufzeiten aufweisen. 

Um die unterschiedlichen Laufzeiten auszugleichen und der Ba- 
sisstation rahmensynchrone Daten liefern zu konnen, misst die 
Mobilstation die Signallauf zeit zur Basisstation und korri- 
giert dadurch einen Sendebeginn ihrer bursts. 

Die Signallaufzeit wird in einem sogenannten „timing advance* 
(TA) codiert und weist eine Abhangigkeit zur Entfernung zwi- 
schen Mobilstation und gesprachsfuhrenden Basisstation auf . 

FUr das TA stehen 64 Stufen zur Verfugung, die mit Werten 0 
bis 63 (bit-) codiert werden und die Laufzeit reprasentieren. 

Da Positionen von Basisstationen bekannt sind, kann aus einem 
TA bzw. aus der Signallaufzeit auf die Position der Mobilsta- 
tion geschlossen werden. 

Eine Messgenauigkeit bei der Lauf zeitbestimmung betragt eine 
Bitdauer, also in GSM 48/13 us, was einer einfachen Weglange 
von etwa 554m entspricht. 

Aus [3] ist eine Positionsermittlung einer mobilen Kommunika- 
tionseinrichtung in einem UMTS-Netz (= Universal Mobile Tele- 
communication System-Net z ) bekannt . 

Bei dem entsprechenden UMTS- Mobilfunkstandard, auf welchem 
das UMTS-Netz basiert, ist eine Positionsermittlung eines Mo- 
bilfunkgerates bereits explizit im Standard enthalten bzw. 
wird durch diesen gefordert (TS 25.305 V3 . 1 . 0 : stage 2 *Func- 
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tional Specification of Location Services in UTRAN" (release 
99), 3GPP TSG-RAN-WG2, 2000). 

Aus [4], [5] und [6] sind weitere Verfahren zur Lokalisierung 
einer mobilen Kommunikationseinrichtung in einem Kommunikati- 
onsnetz bekannt. 

Aus [10] ist ein nichtlineares Filter, ein Prior Density 
Splitting Mixture Estimator (PDSME) , bekannt. 

Dieses nichtlineare Filter, dieser PDSME, basiert auf neue 
Gaussche Mischf ilter-Algorithmen fur nichtlineare Systemzu- 
standsgleichungen und/oder nichtlineare Messgleichungen. 

5 Ublicherweise treten Linearisierungsf ehler auf, wenn die 

nichtlinearen Systemzustandsgleichungen und/oder die nichtli- 
nearen Messgleichungen durch Linearisierungen ersetzt werden 
bei Anwendung von linearen Filtertechniken. 

0 Der PDSME bewirkt eine Reduktion dieser Linearisierungsf ehler 
durch Splitten der a priori Dichte in mehrere Gaussche Misch- 
komponenten mit jeweils kleineren Kovarianzen. 

Der PDSME kann angewandt werden sowohl fur einen sogenannten 
25 prediction Step' bei den nichtlinearen Systemzustandsglei- 
chungen als auch fiir einen sogenannten „Filter Step' bei den 
nichtlinearen Messgleichungen. 

Weiter ist aus [10] bekannt eine Messung fur den Linear isie- 
30 rungsfehler, die vergleichbar einer Kullback-Leibler-Distanz 
ist. Daruber kann eine gewunschte Abschatzungsgtite bei dem 
PDSME eingestellt werden. Durch Beschrankung einer maximalen 
Anzahl von Mischkomponenten bei dem PDMSE kann eine obere 
Grenze fur einen Rechenaufwand angegeben werden. 




<jt\ uj a U3^4y DE 



Die genannten Lokalisierungsverf ahren weisen unter anderem 
die Nachteile auf, dass die durch sie ermittelten Positioner* 
der mobilen Kommunikationseinrichtungen ungenau und demzufol- 
ge mit grofien Unsicherheiten behaftet sind. Genauere Verfah- 
ren erfordern allerdings aufwendige Zusatzeinrichtungen und 
kostspielige Modif ikationen an den bzw. der Kommunikat ions- 
net ze und Kommunikationseinrichtungen. 

Somit liegt der Erfindung die Aufgabe zugrunde, eine genaue 
und mit meglichst geringen Unsicherheiten behaftete Lokali- 
sierung einer mobilen Kommunikationseinrichtung in einem Kom- 
munikationsnetz zu erm6glichen> welche mdglichst einfach und 
kostengiinstig realisiert werden kann. 

Diese Aufgaben werden durch das Verfahren und die Anordnung 
sowie durch das Compute rprogramm mit Programmcode-Mitteln und 
das Computerprogramm-Produkt zur Ermittlung einer ausgewahl- 
ten Position einer mobilen Kommunikationseinrichtung in einem 
Kommunikationsnetz mit den Merkmalen gemafi dem jeweiligen un- 
abhangigen Patentanspruch geldst. 

Bei dem Verfahren zur Ermittlung einer ausgewahlten Position 
einer mobilen Kommunikationseinrichtung in einem Kommunikati- 
onsnetz mit mindestens einer ersten Basisstation, eingerich- 
tet zu einer ersten Kommunikation mit der mobilen Kommunika- 
tionseinrichtung mittels erster Kommunikationssignale, wird 
ein nichtlineares Kommunikationsmodell mit einem determinis- 
tischen Anteil und einem stochastischen Anteil verwendet. 

Dieses nichtlineare Kommunikationsmodell umfasst einen deter- 
ministischen und einen stochastischen Anteil 

- welcher deterministische Anteil eine Abhangigkeit zwischen 
Kommunikationssignale einer Basisstation und Positionen 
der mobilen Kommunikationseinrichtung beschreibt und 
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- welcher stochastische Anteil eine Unsicherheit des deter- 
ministischen Anteils beschreibt. 

Unter Verwendung des zu der ausgewahlten Position zugehorigen 
ersten Konauunikations signals der ersten Basisstation und dem 
nichtlinearen Kommunikationsmodell wird die ausgewahlte Posi- 
tion ermittelt, wobei bei der Ermittlung eine nichtlineare 
Bayesianische Filtertechnik angewandt wird. 

Durch Anwendung der nichtlinearen Bayesianische Filtertechnik 
auf das Kommunikationsmodell bzw. aufi das Kommunikationssig- 
nal wird ein mogliches Auf enthaltsgebiet der ausgewahlten Po ( 
sition bestimmt, aus welchem schlieiilich die ausgewahlte Po- 
sition bestimmt wird. 

Unter der nichtlinearen Bayesianischen Filtertechnik bei der 
Erfindung ist im allgemeinen die Technik wie folgt zu verste- 
hen: 

- die Bayesianische Filtertechnik entspricht einer zeitli- 
chen Fortschreibung bzw. Veranderung (auch Berechnung der 
zeitlichen Veranderung) bekannter Wahrscheinlichkeitsver- 
teilungen der Systemzustande und Systemausgange durch das 
Vorliegen neuer Messwerte. 

Dem erfinderischen Lokalisierungsverf ahren liegt die Idee 
zugrunde, aus in Kommunikationsnetzen, wie zum Beispiel WLAN, 
GSM, DECT-Netzen, zur Verfugung stehenden Kommunikationssig- 
nalen zwischen mindestens einer Basisstation (Sender) und ei- 
ner Mobilstation (Empfanger) entf ernungsrelevante Parameter 
und daraus geographische Inf ormationen, in diesem Fall exn 
mogliches Aufenthalts- oder Entf ernungsgebiet bzw. Unsicher- 
heitsgebiet der Mobilstation, zu gewinnen. 

Anzumerken ist, dass Sender und Empfanger der Kommunikations- 
5 signale technisch auch vertauscht realisiert sein konnen. Die 
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erf inderische Vorgehensweise bleibt davon unbertihrt . In die- 
sem Fall ware die Basisstation der Empfanger und die Mobil- 
station der Sender der in Rede stehenden Kommunikationssigna- 
le. 

5 

Weiter liegt der dem erf inderischen Lokalisierungsverf ahren 
die Erkenntnis zugrunde, dass die von der Mobilstation emp- 
fangenen Kommunikationssignale abhangig sind von der geogra- 
phischen Position der Mobilstation im Kommunikationsnetz bzw. 
10 zu der Basisstation. 

Diese Abhangigkeit zwischen den empfangenen Kommunikations- 
signalen bzw. zwischen Signalstarken der (von der Basisstati- 
on)- empfangenen Kommunikationssignale und der geographischen 
Position der Mobilstation kann durch sogenannte Messgleichun- 
gen, d.h. durch das Kommunikationsmodell, beschrieben werden. 

Diese Messgleichungen bzw. das Kommunikationsmodell umfassen 
einen deterministischen Anteil, welcher deterministische An- 
teil eine Abhangigkeit zwischen dem empfangenen Kommunikati- 
onssignal und der Position der Mobilstation beschreibt, und 
einen stochastischen Anteil, welcher stochastische Anteil ei- 
ne Unsicherheit des deterministischen Anteils beschreibt. 

Diese Unsicherheit kann beispielsweise eine solche des Kommu- 
nikationssignals und/oder eine solche der erwahnten Abhangig- 
keit (Modellunsicherheit) sein . 

Die erf indungsgemafi auf diese Messgleichungen bzw. das Kommu- 
30 nikationsmodell angewandte nichtlineare Bayesianische Filter- 
technik zur Bestimmung des moglichen Auf enthaltsgebiets der 
Mobilstation bzw. der ausgewShlten Position verbessert die 
Qualitat der Lokalisierung verglichen mit Ublichen Filter- 
techniken wie einera Extended Kalmann Filter (EKF) erheblich. 
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Die Anwendung der nichtlinearen Bayesianischen Filtertechnik 
kann in einem sogenannten ^Filter Step" 4 einmalig oder itera- 
tiv mehrmalig erfolgen. 

In diesem moglichen Auf enthaltsgebiet bzw. Unsicherheitsge- 
biet, beschrieben unter Verwendung stochastischer Grofien, wie 
einem Mittelwert und/oder einer Varianz, wird schliefilich die 
Mobilstation bzw. deren Position vermutet . 

Diese Vermutung kann mathematisch formuliert werden durch ei- 
nen Kennwert des moglichen Auf enthalts- bzw. Unsicherheitsge- 
biets, wie einem Schwerpunkt oder einem Erwartungswert . Die- 
ser kann dann als Schatzung fUr die ausgewahlte Position der 
Mobilstation verwendet wird. 

Ein besonderer Vorteil der Erfindung liegt darin, dass die 
Lokalisierung auf Basis von Kommunikationssignalen und be- 
kannten Positionen von Basisstationen durchgeftihrt wird, wel- 
che in einem normalen Betrieb bei einem Mobilfunksystem an- 
fallen und dort zur Verftigung stehen. Dadurch kann auf auf- 
wendige Anderungen und Erweiterungen sowie Zusatzmessungen 
bestehender Mobil funks yst erne bzw. bei bestehenden Mobilfunk- 
systemen verzichtet werden. 

Das Computerprogramm mit Programmcode-Mitteln ist eingerich- 
tet, urn alle Schritte gemafi dem erf indungsgemaJien Verfahren 
zur Ermittlung einer Position, d.h. dem erf indungsgemaBen Lo 
kalisierungsverfahren, durchzufuhren, wenn das Programm auf 
einem Computer ausgeftthrt wird. 

Das Computerprogramm-Produkt mit auf einem maschinenlesbaren 
Trager gespeicherten Programmcode-Mitteln ist eingerichtet, 
urn alle Schritte gemafi dem erf indungs gemafi en Lokalisierungs- 
verfahren durchzufuhren, wenn das Programm auf einem Computer 
ausgeftlhrt wird. 
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Die Anordnung sowie das Computerprogramm mit Programmcode- 
Mitteln, eingerichtet urn alle Schritte gemaB dem erfinderi- 
schen Lokalisierungsverf ahren durchzuftihren, wenn das Pr,o- 
gramm auf einem Computer ausgefiihrt wird, sowie das Computer- 
programm-Produkt mit auf einem maschinenlesbaren Trager ge- 
speicherten Programmcode-Mitteln, eingerichtet um alle 
Schritte gemaB dem erf inderischen Lokalisierungsverf ahren 
durchzuftihren, wenn das Programm auf einem Computer ausge- 
fiihrt wird, sind insbesondere geeignet zur Durchftihrung des 
erf indungs gemaB en Lokalisierungsverf ahrens oder einer seiner 
•nachfolgend erlauterten Weiterbildungen. 

Bevorzugte Weiterbildungen der Erfindung ergeben sich aus den 
abhangigen Anspruchen. 

Die im weiteren beschriebenen Weiterbildungen beziehen sich 
sowohl auf die Verf ahren als auch auf die Anordnung. 

j 

Die Erfindung und die im weiteren beschriebenen Weiterbildun- 
gen konnen sowohl in Software als auch in Hardware, bei- 
spielsweise unter Verwendung einer speziellen elektrischen 
Schaltung, realisiert werden. 

Ferner ist eine Realisierung der Erfindung oder einer im wei- 
teren beschriebenen Weiterbildung moglich durch ein computer- 
lesbares Speichermedium, auf welchem das Computerprogramm mit 
Programmcode-Mitteln gespeichert ist, welches die Erfindung 
oder Weiterbildung ausfuhrt. 

Auch kann die Erfindung oder j ede im weiteren beschriebene 
Weiterbildung durch ein Computerprogrammerzeugnis realisiert 
sein, welches ein Speichermedium aufweist, auf welchem das 
Computerprogramm mit Programmcode-Mitteln gespeichert ist, 
welches die Erfindung oder Weiterbildung ausfUhrt. 
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Bei einer Kommunikation in einem Kommunikationsnetz zwischen 
einer mobilen Kommunikationseinrichtung (Mobilstation) , bei- 
spielsweise einem Mobiltelef on, und einer Basisstation, bei- 
5 spielsweise einer Rundantenne bzw. einem Rundstrahler Oder 
einer bzw. mehrerer sektoraler Antennen, werden Daten, die 
Kommunikationssignale, in Signalpaketen, sogenannten bursts, 
tibertragen. 

10 Basierend auf den bzw. unter Verwendung der tibertragenen Kom- 
munikationssignalen bzw. Signalpaketen lassen sich verschie- 
dene entf ernungsrelevante Parameter ermitteln, welche wieder^ 
um als Grundlage fur die Ermittlung von den moglichen Aufent* 
halts- bzw. Entf ernungsbebieten herangezogen werden konnen. 

15 

Ein solcher entf ernungsrelevanter, d.h. entf ernungsabhangi- 
ger, Parameter ist beispielsweise eine Feldstarke eines Sig- 
nalpakets . 



20 Die Feldstarke weist eine nattlrliche Abhangigkeit zur Entfer- 
nung zwischen der Mobilstation und der (gesprachsftihrenden) 
Basisstation auf und liefert demzufolge eine Information tiber 
das mogliche Auf enthaltsgebiet bzw. Entf ernungsgebiet (Unsi- 
cherheitsgebiet) der Mobilstation. 



25 




Diese Abhangigkeit zwischen Feldstarke und Entfernung kann 
durch physikalische Modelle, welche ein Ausbreitungsverhalten 
von Signalen beschreiben, beschrieben werden. 



30 Ftir die Ermittlung der ausgewahlten Position kann ein Kenn- 

wert des m6glichen Aufenthalts- bzw. Unsicherheitsgebiets be- 
stimmt werden, wie ein Schwerpunkt oder ein Erwartungswert, 
welcher dann als Schatzung ftir die ausgewahlte Position der 
Mobilstation verwendet wird. 



35 
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Die Qualitat der Lokalisierung kann weiter verbessert werden 
unter Verwendung eines Benutzermodells, welches eine Bewegung 
der mobilen Kommunikationseinrichtung beschreibt. So kann ein 
solches Benutzermodell eine maximale Schrittweite bzw. Bewe- 
gungsdistanz fUr einen vorgegebenen Zeitschritt begrenzen. 

Die Erfindung eignet sich insbesondere zu einem Einsatz im 
Umfeld eines digitalen, zellularen Mobil funksyst ems , wie ei- 
nes GSM-Netzes, und dort beispielsweise zur Lokalisierung ei- 
nes GSM— Telefons (Mobiltelef on) . 

Dabei werden bei dem Einsatz der Erfindung nur die dem Mobil- 
telefon zur Verftigung stehenden Daten verwenden, wobei weder 
am GSM-Netz noch an Mobilstationen in dem GSM-Netz kostspie- 
lige Anderungen vorzunehmen sind. 

Beispielsweise sind bekannt von einem GSM-Netz die Positionen 
der einzelnen Basisstationen und deren Antennen sowie derei* 
Charakteristika, welche Auskunft uber das Versorgungsgebiet . 
der jeweiligen Antenne gibt. 

Das Mobiltelefon seinerseits steht ftir einen korrekten Ver- 
bindungsaufbau immer in Kontakt mit den 

empfangbaren Antennen, urn die am besten far ein Gesprach 
geeignete Antenne vom Netz zugeteilt zu bekommen. Dazu misst 
es unter anderem die Empf angsf eldstarke der empfangbaren An- 
tennen sowie bestimmt Signallauf zeiten, die dann ebenfalls 
bekannt sind, 

Auf Basis dieser vorliegenden Inf ormationen erfolgt dann die 
Lokalisierung des Mobiltelef ons gemafi der erf inderischen Vor- 
gehensweise . 

Auch eignet sich die Erfindung zu einem Einsatz im Umfeld 
weiterer digitaler, zellularen Mobilfunksysteme, wie eines 
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WLAN Oder eines DECT-Netzes [11], und dort beispielsweise zur 
Lokalisierung eines DECTHYIobiltelef ons • 

In Figuren ist ein Ausf Uhrungsbeispiel der Erfindung darge- 
5 stellt, welches im weiteren naher erlautert wird. Das Ausf uh- 
rungsbeispiel ist dabei unterteilt in einen grundlegenden 
Teil, welcher Grundlagen der erf inderischen Vorgehensweise 
darstellt, und einen angewandten Teil/ welcher anhand eines 
konkreten Zahlenbeispiels die erf inderische Vorgehensweise 
10 sowie Ergebnisse verdeutlicht bzw. angibt. 

Es zeigen 

Figur 1 eine Skizze, in welcher eine angenommene lineare Na- 
15 herung fur das nicht lineare Kommunikationsmodell dar- 

gestellt ist; 

Figur 2 eine Skizze, in welcher die Unsicherheit der Kommuni- 
kationssignale durch rauschen dargestellt ist; 

20 

Figur 3 eine Skizze, in welcher ein Oberblick tiber den PDSME- 
Algorithmus mit einem linearisierten Fiterstep (obe- 
rer Ausschnitt) und einem Prognoseschritt (Predikti- 
onsstep) (unterer Ausschnitt) dargestellt ist; 

25 

Figuren 4a und b Skizzen, in welchen ein Evaluierungs- 

Szenario ftir die Lokalisierung gemaii dem AusfUhrungs- 
beispiel dargestellt is; 

30 Figuren 5a und b Skizzen, in welchen eine Lokalisierung unter 
Verwendung eines Extended Kalman Filters dargestellt 
ist (erster Filterschritt Fig. 5a; 170-ter Filter- 
schritt Fig, 5b) ; 
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Figuren 6a und b Skizzen, in welchen eine Lokalisierung unter 
Verwendung des PDSME gemaii dem Ausf uhrungsbeispiel 
dargestellt ist (erster Filterschritt Fig. 6a; 170- 
ter Filterschritt Fig. 6b; 

Figuren 7a und b Skizzen, in welchen Ergebnisse der EKF- 

Filterung und der PDSME- Fi It erung gemafi dem Ausf uh- 
rungsbeispiel in Bezug auf die korrekten Ergebnisse 
dargestellt sind; 

Figur 8 eine Skizze, in welcher eine Anzahl von Gausschen 
Mischkomponenten in einer Abhangigkeit von Filter- 
schritten dargestellt sind. 

Ausftihrungsbei spiel : Lokalisierung eines DECT-Mobiltelef ons 
in einem DECT-Netz basierend auf einer nichtlinearen Filter- 
technik 

1 . Grundlagen 

Nachfolgend wird die Anwendung einer nichtlinearen Bayesiani- 
schen Filtertechnik bei der Lokalisierung eines mobilen Funk 
telefons beschrieben. Die Anwendung wird beschrieben am Bei- 
spiel einer Lokalisierung eines DECT-Mobiltelef ons in einem 
DECT-Netz mit mehreren Basisstationen. 

Signals tarken (Feldstarken) der von einem Mobil telef on emp- 
fangenen und dort gemessenen Kommunikationssignale (einer 
diese Signale aussendenden Basisstation) sind abhangig von 
der Position des Mobiltelef ons bezttglich der die Kommunikati 
onssignale aussenden Basisstation. Dieser Zusammenhang wird 
beschrieben durch nicht lineare Messgleichungen (nichtlinea- 
res Kommunikationsmodell) . 
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Diese nicht linearen Messgleichungen bzw. dieses nicht linea- 
re Kommunikationsmodell umfasst einen deterministischen An- 
teil, welcher die empf angenen Signals tar ken als Funktion der 
Position beschreibt, und einen stochastischen Anteil, welcher 
Modellfehler und Messrauschen berttcksichtigt . 



Zusatzlich werden Benutzermodelle mit berticksichtigt, welche 
ein Wissen tiber eine (raumliche) Bewegung von Benutzern der 
Mobiltelefone mit einbringen. 



Eine neue nichtlineare Filtertechnik, ein sogenannter „Prior 
Density Splitting Mixture Estimator* (PDSME) , einordnbar als 
Gausscher Mischf ilteralgorithmus, verbessert signifikant eid 
Lokalisierungsqualitat verglichen mit Standardf iltertechni- 
ken, wie ein ^Extended Kalman Filter" (EKF) (vgl. Figuren 5 
bis 7) . 




Typische Anwendungen ftir Lokalisierungsauf gaben sind Positi- 
onsermittlungen mobiler Kommunikationseinrichtungen (Mobilge- 
raten) in WLANS/ GSM-Netzen oder DECT-Netzen. 

FUr die Positionsermittlung sollten dabei Inf ormationen ver- 
wendbar sein, welche bereits wahrend einer normalen Operati- 
onsphase, d.h. bei einer „normalen* Kommunikation zwischen 
Mobilgerat und Basisstation, zur Verftigung stehen. I 

Das aus [12] bekannte Lokalisierungsverf ahren fUr Mobiltele- 
fone in GSM-Netzen basiert auf ortsabhangige, von den Mobil- 
telefonen empfang- und dort messbare Signalst&rken der von 
Basisstationen ausgesendeten Kommunikationssignale . Wegen der 
Ortsabhangigkeit sind diese charakteristisch fttr die Position 
des Signalempfangers, in diesem Fall fUr die Position des Mo- 
biltelefons bzw. eines Benutzers des Mobiltelef ons . 

Wahrend der normalen Operationsphase werden die empfangenen 
Signalstarken aller empfangbarer Kommunikationssignale der 
jeweiligen Sender (Basisstationen) durch das Mobiltelef on ge- 
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messen und verglichen. Dadurch wird ein sogenanntes Handover 
des Mobiltelefons zwischen verschiedenen Sendern bzw. Basis- 
stationen gesteuert . 

Basierend auf diese im Normalbetrieb anfallenden Messungen 
konnen Lokalisierungsverf ahren entwickelt werden, urn zusatz- 
lich die Position des Mobiltelefons in Bezug auf die Sender 
bzw. im GSM-Netz zu ermitteln. 

Aus [13] und [17] sind weitere Vorgehensweisen zur Lokalisie- 
rung mobiler Kommunikationseinrichtungen bekannt. Diese Vor- 
gehensweise bei der Lokalisierung basieren auf einer Anwen- 
dungen von nichtlinearen Zustandsschatzungen. 

Nachfolgend wird die Anwendung einer neuen nichtlinearen Fil- 
tertechnik fttr die Lokalisierung von DECT-Mobiltelef onen be- . 
schrieben. 

Um die Lokalisierung von Mobiltelef onen zu ermoglichen, wird, 
wie aus [18] bekannt ist, ein stochastisches Naherungsmodell * 
far die Ausbreitung von Radiowellen Oder empf angbaren Signal- 
starken von Kommunikat ions signal en in einer Lokalisierungsum-: 
gebung identif iziert bzw. entwickelt. 

Ein entsprechendes Ausbreitungsmodell kann auch grundsatzlich 
durch physikalische Ausbreitungsmodelle von elektromagneti- 
schen Wellen erlangt werden. 

Eine Anpassung rein physikalischer Ausbreitungsmodelle an re- 
ale (Lokalisierungs-)Umgebungen ist jedoch sehr kompliziert, 
da Reflexionen oder Interf erenzen oder andere physikalische 
Phanomene die physikalisch angenommene Ausbreitung beeinf lus- 
sen \ind verf^lschen. 

Da reale Ausbreitungsparameter , welche die elektromagneti- 
schen (Ausbreitung-) Eigenschaf ten in der realen Lokalisie- 
rungsumgebung, in welcher das Mobiltelefon zu lokalisieren 
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ist, beschreiben, nur partiell bekannt bzw. bestimmbar sind, 
ist die Anwendung solcher physikalisch basierter Ausbrei- 
tungsmodelle fttr solche realen Anwendungen ausgeschlossen. 

5 Deshalb ist ein Ausbreitungsmodell, welches auf tatsachlichen 
Messungen der (logarithmierten) empfangenen Signals tarken der 
Kommunikationssignale in der tatsachlichen Lokalisierungsum- 
gebung beruht bzw. sich daraus ableiten lasst, theoretischen 
physikalischen Ausbreitungsmodellen vorzuziehen. 

10 

Die Entwicklung eines solchen aus tatsachlichen Messungen ab- 
geleitetes „Messmodell* bzw, Kommunikation- bzw. Ausbrei- m 
tungsmodell erfordert ein Einmessen oder eine Kalibrierung irf 
15 Rahmen dessen Parameter des Messmodells festgelegt werden. 

Einmessen/Kalibrieren („model generation phase*) 

Ftlr ein solches Einmessen/Kalibrieren wird der Empfanger (be- 
20 nutzt im Sinne eines Messgerates) , d.h. das Mobiltelef on, an 
mehrere Positionen in der Lokalisierungsumgebung bewegt, bei- 
spielsweise an Gitterpunkten erhaltlich durch Platzierung ei- 
nes Gitters auf die Lokalisierungsumgebung. 

25 An jedem Gitterpunkt werden die Feldstarken aller dort emp- 
fangbaren Basisstationen gemessen und in einer Karte der Lo- 
kalisierungsumgebung gespeichert. Auf diese Weise wird die 
Feldstarkenverteilung jeder Basisstation ttber die Lokalisie- 
rungsumgebung bestimmt. Die Feldstarkenverteilungen konnen 

30 durch das Messmodell modelliert werden. 

Basierend auf diese Modellinf ormationen (und von an einer ak- 
tuellen Position gemessenen Feldstarken) kann dann die aktu- 
elle, von einem Empfanger bzw. von dem Mobiltelef on in der 
35 Lokalisierungsumgebung eingenommenen Position bestimmt werden 
(Lokalisierung in einer ^localisation phase* ) . 
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Weiter wird nachfolgend beschrieben werden, wie ein sto- 
chastisches Messmodell in der model generation phase identi- 
fizierbar ist. 

Dieses stochastische Messmodell besteht aus einer analyti- 
schen, deterministischeh Messf unktion, welche die logarith- 
mierte, empfangene Signalf eldstarke in Abhangigkeit von Posi- 
tionskoordinaten beschreibt. Dartiber hinaus wird ein sto- 
chastischer Anteil (Komponente) identif iziert , welcher Mo- 
dellunsicherheit wie auch Messrauschen berUcksichtigt. 

Ein solcher stochastischer Ansatz bei der Lokalisierung be- 
rUcksichtigt die Unsicherheiten bei der ermittelten Position 
durch Wahrscheinlichkeiten bzw. Wahrscheinlichkeitsdichte- 
funktionen. 

Neben dem Messmodell ftir die empf angbaren Kommunikationssig- 
nale wird die GUte der Lokalisierung dadurch weiter verbes- ; 
sert, dass ein Benutzermodell mit eingebracht wird. Ober die- 
ses Benutzermodell kann Wissen Ober eine (raumliche und/oder 
zeitliche) Bewegung eines Benutzers mit berUcksichtigt wer- 
den. 

Durch solche Benutzermodelle konnen maximale „Schrittweiten* 
von Benutzern in einem Zeitschritt, welcher Uber zwei zugeho- 
rige bzw. korrespondierende Signalmessungen (am Anfang und am 
Ende des Zeitschritts) definiert wird, beschrankt ist. 

Ein zusatzliche Verbesserung der Lokalisierungsgute lasst 
sich erreichen durch rekursive Positionsermittlung durch Kom- 
bination mehrfacher verschiedener nacheinanderf olgender Mes- 
sungen . 
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Der hier benutzte neue stochastische Ansatz bei der Lokali- 
sierung fuhrt zu nichtlinearen, mehrdimensionalen Messglei- 
chungen. Eine exakte Losung dieses Bayesianischen Filterprob- 
lems fiihrt zu komplizierten nicht-Gausschen Wahrscheinlich- 
keitsdichten, welche die gesuchte Position mit ihren Unsi- 
cherheiten beschreiben. 



10 



25 



Diese Wahrscheinlichkeitsdichte (-funktionen) kOnnen durch 
Gaussche Mischdichten [19], [20] angenahert werden. 



Diese Annaherung wird hier durch einen neuentwickelten soge- 
nannten „Prior Density Splitting Mixture Estimator* (PDSME) 
durchgefUhrt, der entwickelt wurde, urn eine durch einen An- 
wender einstellbare Gute der Lokalisierung zu ermoglichen - 
15 im Gegensatz zu tiblichen, nichtlinearen, bekannten Filter- 
techniken, wie einem Extended Kalman Filter (EKF) . 

Jedoch kann eine Angabe einer oberen Grenze ftir einen Rechen- 
aufwand durch BeschrSnkung der Anzahl von anwendbaren Gauss- 
20 chen Mischkomponenten die Lokalisierungsgute reduzieren. Dar- 
tlber hinaus wird durch die Representation von Dichten bei den 
Gausschen Mischungen eine rekursive Vorgehensweise sowohl bei 
den nichtlinearen ..Predictions Steps* als auch bei den nicht- 
linearen „Filtersteps* moglich. 





Weiter wird im ersten Teil beschrieben werden, wie das sto- 
chastische Lokalisierungsproblem formuliert werden kann. Die 
deterministischen Anteile und die stochastischen Anteile des 
Messmodells werden nachfolgend erklart werden. Die Anwendung 
30 von Bayesianischen Filteralgorithmen, basierend auf eine 

gausschen MischnSherung ftir die exakten Wahrscheinlichkeits- 
dichten verwendet in dem neuen PDSME wird beschrieben werden. 



35 



Im zweiten Teil werden Ergebni sse einer Lokalisierung eines 
DECT-Mobiltelefons in einem DECT-Netz dargestellt und die Gti- 
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te der Lokalisierung unter Verwendung des PDMSE mit der unter 
Verwendung eines bekannten . EKF verglichen werden. 

Formulierung des Lokalisierungsproblems 

Die Lokalisierungsaufgabe fur ein mobile Kommunikationsein- 
richtung in einem Kommunikationsnetz kann in zwei Grundprob- 
leme eingeteilt werden: 

1) die Identifikation des Messmodells mit einem deterministi- 
schen und einen stochastischen Anteil in der model genera-, 
tion phase; 

2) die Ermitt lung. der Wahrscheinlichkeitsdichte der gesuchten 
Position in einem Bayesianischen Filterstep in der Lokali- 
sierungsphase . 

Das Messmodell 



beschreibt die Messung der logarithmierten empfangenen Sig- 
nalstarken 



des C-then Transmitters (Senders), Q= 1, . . . , N, als eine de- 
terministische nichtlineare Funktion hC(xk) der Positionsko- 
ordinaten xk=[xl / k / x2,k]far jeden der N verschiedenen Sen- 
der. 

Der stochastische Anteil vk des Messmodells wird reprasen- 
tiert durch eine additive Unsicherheit . Additive Unsicherhei- 
ten im Messmodell der logarithmierten empfangenen Signalstar- 
ke entsprechen multiplikativen Unsicherheiten der Signalleis- 
tung P, welche auf einen Einfluss von unbestimmten Dampfungen 
zuruckzufuhren sind. 




(1) 
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Jeder Transmitter kann identif iziert werden Uber eine ihm ei- 
genen ID, welche durch die Basisstation wahrend einer Kommu- 
nikation tibertragen wird. 



Dadurch wird das Lokalisationsproblem erheblich vereinfacht, 
weil jeder Messwert direkt der zugehorigen nichtlinearen 
Messgleichung zugewiesen werden kann. 

In der Lokalisierungsphase werden die nichtlinearen Filteral- 
gorithmen in jedem Zeitschritt k angewandt. Jedoch k6nnen oh- 
ne die additiven Unsicherheiten vk die impliziten, nichtline- 
aren Beziehungen zwischen dem Vektor mit den gemessenen loga 
rithmierten empfangenen Signalstarken yk der Transmitter una 
den Positionen xk des Empfangers nur numerisch gelost werden, 
beispielweise durch eine Methode der kleinsten Fehlerquadra- 
te. Berticksichtigt man die Unsicherheiten mit, so konnen die 
Messgleichungen verwendet werden, urn die gesuchte Position im 
Filterstep zu erneuern („update*). 

Die weiter beschriebene Vorgehensweise bei der Lokalisierung 
basiert auf einer Naherung der exakten Losung des bayesiani- 
schen Filtersteps durch Gaussche Mischdichten 




= 2>& — * 1 ' j ' 



mit L Komponenten, definiert durch nicht-negative Gewichte w, 
Erwartungswerte |i und Kovarianzen C. 

1st der Messwert y verfUgbar, so ergibt sich die exakte Lo- 
sung des Bayesianischne Filtersteps aus : 

fjfepkr C) = c f x p (x k ) fv, ^(Yk, ? - h^(x k )) 

mit der nichtlinearen Messfunktion h(x), den additiven Unsi- 
cherheitsdichtefunktiuonen f (v) und eine Normalisierungskon- 
stante c. 
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A priori Wissen Ube die gesuchte Position wird repr^sentiert 
durch die Dichtefunktion f(x). 

Fur jede neue Messung wird die vorangehende a posteriori 
Dichte f (x|y) interpretiert als neue a priori Dichte, so dass 
die gesuchte bzw. ermittelte Position rekursiv erneuert (up- 
date) werden kann. 

Analog zum Bayesianischen Filterstep wird ein Predikti- 
onsstep, welcher zur Beschreibung des nichtlinearen Benutzer- 
modells 

Sc+i = a k (£k) + (2) 

verwendet wird, ebenso zur Naherung der exakten Dichte durch 
Gaussche Mischungen durchgef tihrt . 

Die nichtlineare Funktion a(x) ist ein deterministisches Mo-, 
dell einer Benutzerbewegung. Unsicherheiten werden wieder 
durch additives Rauschen w beriicksichtigt . 

Stochastische Modellierung 

Das Messmodell wwird identif iziert tiber Messungen der loga- 
rithmierten Signalstarken ftlr jeden Transmitter an Gitter- 
punkten eines Gitters,, welchen die Lokalisierungsumgebung 
abdeckt. Die deterrainistische (h(x))und die stochastische 
Komponente (v) des Messmodells (1) werden vor der Lokalisie- 
rung identif iziert . 

Der deterministisch, analytische Anteil h(x) ist eine Mess- 
gleichung, welche die empfangene Signalleistung als eine 
Funktion der Positionskoordinaten beschreibt. Der stochasti- 
sche Anteil v ist ein Modell fur die Unsicherheiten des de- 
terministischen Antteils. Diese Unsicherheiten setzen sich 
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zusammen aus raumlichen Unsicherheiten entsprechend Modell- 
fehlern und zeitliches Messrauschen. 

Deterministisches Messmodell der logarithmierten empfangenen 
5 Signalstarken 

Fig.l zeigt, dass die Abnahme der logarithmierten, empfange- 
nen Signalstarke Uber eine Distanz von mehreren Metern anna- 
hernd linear ist. In anderen Dokumenten werden ahnliche An- 
10 nahmen beziiglich des Messmodells als „linear-loss-model* [24] 
Oder „linear-scope-model* [18] bezeichnet. 

Mathematisch lasst sich fur zweidimensionale Positionskoordi 
naten x richtungsabhangige lineare Abnahme der logarithmier- 
15 ten empfangenen Signalstarken beschreiben durch N unabhangige 
Messgleichungen 

Die Parameter m£ und PC der positiven semidef initen quadrati- 
schen Form und der additive Offset AC werden fur jeden Trans- 
20 mitter C= 1/ . N identif iziert . 

Die 6N Parameter mC, PC und AC werden bestimmt durch eine Me- 
thode der kleinsten Fehlerquadrate, bei welcher die Abwei- 
chung zwischen dem angenaherten Messmodell und der gemessenen 
25 logarithmierten Signalstarken auf den Gitterpunkte, gemessen 
wahrend der model generation phase , minimiert wird. Zur Re- 
duktion des zeitlichen Messrauschens wird der Mittelwert tiber 
mehrere Messungen an jedem Gitterpunkt zur Berechnung des 
Durchschnitts der gemessenen Werte verwendet. 

30 

Stochastische Modellierung der Unsicherheiten 



Das stochastische Unsicherheitsmodell berucksichtigt sowohl 
die Abweichung zwischen den angenaherten, im obigen beschrie- 
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benen Modell und der wahren Verteilung der logarithmierten 
empfangenen Signalstarken, gemessen auf den Gitterpunkten, 
und das zeitliche Messrauschen. 

5 Modellunsicherheiten 

Bei der beschriebenen Vorgehensweise wird die Abweichung zwi- 
schen den wahren logarithmierten empfangenen Signalstarken 
und dem deterministischen Messmodell h(x) beschrieben durch 
Gaussche Rauschdichten mit entsprechenden Mittelwerte und 
Standardabweichungen fUr jeden Transmitter. Diese Gaussche 
Unsicherheit ist eine Reprasentation far einen Naherungsf eh- 
ler des deterministischen Anteils des Messmodells tiber die 
gesamte Lokalisierungsiamgebung . 
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Time Series 1 
Time Series 2 
Time Series 3 



fgure 2. Spatial variations of the logarithmic received signal power over short distances and temporal measurement 
noise (three different time series). 

3.2.1. Model uncertainties 

In this localization approach, the deviation between the true logarithmic received signal power and the determin- 
istic measurement model h (x k ) is described by a Gaussian noise density with mean standard deviation 
°"^C for eax ^ transmitter. This Gaussian uncertainty is a representation for the approximation error of the de- 
terministic component of the measurement model over the whole localization environment. First, it represents 
the mean deviation over the localization environment between the model and the true received power, because of 
incorrect assumptions for the deterrninistic part of the measurement equation. The uncertainties resulting from 
large-scale deviations between the true received signal power and the measurement model are depicted in Fig. 1. 
Second, measurements have also shown that there exist local deviations, that might be caused by reflexion, 
non-homogenous propagation of the radio waves and interference. These spatial variations of the logarithmic 
received signal power can be seen in Fig. 2, where measurements have been collected with a distance of 2 cm. 

3.2.2. Measurement noise 

In addition to spatial variations of the received signal power, temporal measurement noise can be determined by 
analyzing a series of several different measurements at a fixed position. This temporal measurement noise is also 
shown in Fig. 2, for three different measurements of the logarithmic received signal power at each measurement 
position. This noise is again approximated by a Gaussian density with mean ^ and standard deviation a^. 

3.2.3. Combination of model uncertainties and measurement noise 

lb obtain a simple model comprising both uncertainties described in 3.2.1 and 3.2.2, it is further assumed, that 
both uncertainties are independent. Therefore, they can be modeled by a single Gaussian density f Vt ^ (v^) for 
each transmitter £ = 1, . . . , JV, defined by the mean 





and the standard deviation 



a > - iz< 1 > -i- ia< 2 > 



This model not only implies independent uncertainties, but it is also assumed that both uncertainties can be 
described without considering any position dependency. 




4. FILTERING ALGORITHM 

In this Section, a brief overview of the Prior Density Splitting Mixture Estimator (PDSME) used for the lo- 
calization of radio communication devices is given. Additionally, an adaptation of the measurement equations 
hfefc) derived in Section 3 is introduced to simplify the calculation of the PDSME. Furthermore, a prediction 
step for a simple user model is presented. 

4.1. The PDSME measurement update 

The PDSME algorithm for the localization of radio communication devices presented in this paper is based on 
the calculation of a linearized measurement update for Gaussian mixture densities. The measurement update 
step of this filtering algorithm is shown in a block diagram in the upper part of Fig. 3. Splitting is based on the 
calculation of the linearization error 

R 2 

for each component of the posterior Gaussian mixture density. This criterion is very similar to the Kullback- 
Leibler distance 14 




between the exact posterior density /J* (z k \Vk£) a?d its approximation / x " {* k \yk£) by replacing the nonlinear 
measurement equation h(x k ) by its linearization h {x^) at the mean of the i-th component of the prior density 
function. Calculating the linearization error (4) , the prior Gaussian mixture components which contribute most to 
the approximation error of the posterior density are identified, lb reduce this linearization error, these Gaussian 
mixture components are replaced by splitting them into several mixture (components with smaller covariances 
using splitting libraries, which have been optimized off-line (see Fig. 4 in Jf6 ). L*6j 

After this "analytic resampling" of the prior density, the filter step can be calculated by a bank of EKFs, 
linearizing the measurement equation at the mean of each component of the Gaussian mixture representation of 
the prior density. 

Afterwards, a merging step reduces redundancy by combining several Gaussian mixture components to a 
single Gaussian with negligible approximation error. 

4.2. Adaptation of the measurement model for simplifications of the linearization error 
criterion 

For Gaussian measurement noise, the linearization error (4) can be calculated analytically as a linear combina- 
tion of moments of the density J^fefcl^.c) for polynomial measurement equations h c (x k ). Therefore, in this 
Subsection, an adaptation of the measurement equation (3) is derived. After some algebraic conversions, the 
squared measurement equation 

(5) 




can be re-written as a polynomial function with a modified "measured value" z fc ,< and the transformed uncertainty 

Because of the nonlinear transformation of the random variable v c , the probability density function /a.c 
is no longer Gaussian. In the localization experiment in Section 5, the exact first and second order moments 
of V( are computed to determine a Gaussian approximation of /tj,c C^c)- Note, that the moments of vc are 
depending upon the measured value fax- Therefore, they have to be re-computed for each new measurement of 
the logarithmic received signal power and are not time-invariant as the parameters /* v ,< and <?vX described in 
Section 3. 
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Figure 3. Overview of the PDSME algorithm: The PDSME consists of a linearized filter step (upper part) and linearized 
prediction step (lower part). 



4.3. User modeling by PDSME prediction step 

Similarly to the filter step, the PDSME can also be applied to the calculation of nonlinear prediction steps. 
Analogously to nonlinear filter steps, the calculation of an approximated prediction step also consists of the 
evaluation of a linearization error, a bank of linearized prediction steps and the reduction of the number of 
Gaussian mixture components in a merging step (see the lower part of Fig. 3). In this paper, only a linear user 
model is considered. Therefore, the prediction step can be calculated analytically, since the posterior density has 
been approximated by a Gaussian mixture density in the filter step. 

The prediction model consists of the linear state equation 



£jfc+i =£fc +mk i 
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where the mean jx^ of the additive system noise represents knowledge about possible directions and mean 
step lengths of the user's movement. The covariance matrix of 20* specifies an estimate for the distribution of 
the user's step lengths. For each component i = 1, . . . , L of the Gaussian mixture density, the predicted Gaussian 
mixture component is then described by the mean 

and the covariance 

The estimated position can then be calculated as a weighted superposition 

h 




of the means of ail Gaussian mixture components. 

Till 2 ' 5. LOCALIZATION EXPERIMENT 

In this Section, a real-world localization experiment is presented for the validation of the described localization 
approach for DECT mobile telephones. In the Subsections 5.2 and 5.3, the superior performance of the PDSME 
compared to the Extended Kalman Filter (EKF), 16 a widely used standard approach for the state estimation of 
nonlinear systems, is shown. 

5.1. Evaluation scenario 

In this localization experiment, N = 10 transmitters have been placed in an indoor area of approximately 
30 m x 30 m in one floor of a building. In the model generation phase, the logarithmic received signal power 
of each transmitter has been measured on a grid with 1 m distance between grid points. N = 10 measurement 
equations h{ altogether consisting of 60 parameters for the deterministic components and 20 parameters 
for the stochastic components v k , have been identified. Except for a few areas, influenced by high attenuation 
of the radio waves because of ferroconcrete walls, the assumed model is an appropriate approximation of the 
distribution of the received signal power. 

The initial probability density of the position is chosen as a Gaussian density with the initial mean 



and the initial covariance 




i.e., almost no prior knowledge about the position is available. 

In the localization phase, a measurement of each receivable transmitter is collected along a line as shown 
in Fig. 4, with a distance of At = 0.5m between the true measurement positions P lt P 2 > • • • i ^17- The position 
coordinates of the measurements have been determined to compare the ground truth to the estimated positions. 
Neither in the localization approach by the PDSME nor by the EKF, the exact measurement positions have been 
used for the localization of the receiver. After the measurement update has been calculated for each receivable 
transmitter in a fixed position P< , i = 1, . . . , 17, a prediction step according to Subsection 4.3 has been calculated. 
In this example, the simple user model is defined by the mean 

pi = Om 

and the covariance 



-r. 3 



i.e., no knowledge about a preferred direction of the motion of the user is available. Only the distribution of the 
user's step lengths is represented by the covariance matrix C u . 
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Figure 4. Evaluation scenario of the localization experiment, 
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5.2. Extended Kalman Filter (EKF) 

Using the EKF, the measurement equation (5) is lmearized at the mean of the prior density function. In Fig. 5, 
the results of the EKF are presented for the first and the 170th filter step. In each filter step, the measurement 
of the logarithmic received signal power of a single transmitter is used to update the estimated position. 




Figure 5. Localization with the EKF: First filter step (left) and 170th filter step (right). 

It can be noticed, that there is a significant estimation error after the first filter step. The true position 
of the receiver is not within the support of the estimated posterior density function depicted by its contour 
plot. Furthermore, there is no intersection between the true measurement equation and the estimated density. 
Therefore, it is not possible to apply data validation techniques to find out whether a measured value can be 
"explained" by the estimated density. 



5.3. Prior Density Splitting Mixture Estimator (PDSME) 

In Fig. 6, the posterior densities estimated by the PDSME algorithm are also shown for the first and 170th filter 
step. Obviously, the approximation of the non-Gaussian posterior density in the first filter step, which is very 
close to the numerically calculated optimal Bayesian solution of the filter step, is much better than in the case 
of the EKF. Hence, data validation techniques can now be successfully applied. 
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Figure 6. Localization with the PDSME: First filter step (left) and 170th filter step (right). 

In Fig. 7, the expected values of the posterior density functions calculated by the EKF and the PDSME 
are compared to the true position. Comparing this figure to Fig. 8, it can be seen, that for almost; linear 
filtering problems, i.e., if the covariance of the estimated position is small compared to the nonlinearity of the 
measurement equation after several filter steps, the EKF and the PDSME yield almost the same results. In 
these cases, the PDSME only uses a moderate number of Gaussian mixture components, whereas for strong 
nonlinearities at the begirming of the localization experiment a higher number of approximation components is 
necessary to reduce the estimation error. The superior performance of the PDSME compared to the EKF is also 
shown by the average estimation error 

F k=i 

over the N F = 170 filter steps, which is 3.30 m for the EKF and 1.22 m for the PDSME. 

6. CONCLUSIONS 

In this paper, a stochastic approach for the localization of radio communication devices has been presented, 
which is based on measuring the logaritlimic signal power of the receivable transmitters by a mobile part. A 
stochastic measurement model, consisting of a deterministic and a stochastic component has been identified for 
each transmitter. This measurement model has been used for the estimation of the position of the receiver by 
a novel Gaussian mixture estimator, which is based on splitting the prior density according to a linearization 
error criterion. This criterion is very similar to the Kullback-Leibler distance between the true and the approx- 
imated posterior density, calculated by a linearization of the measurement equation. Applying this estimation 
technique to the localization of DECT mobile telephones, significant improvements of the estimation quality 
can be achieved, if the PDSME is used instead of standard approaches like the EKF. Further improvement 
of the localization quality can be achieved by identifying better detenninistic measurement models and more 
precise characterizations of the measurement noise, which do not assume independence between the different 
uncertainties mentioned in this paper. 
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Pa ten tan spriiche 

l.Verfahren zur Ermittlung einer ausgewahlten Position einer 
mobilen Koramunikationseinrichtung in einem Kommunikationsnetz 
5 mit mindestens einer ersten Basisstation, eingerichtet zu ei- 
ner ersten Kommunikation mit der mobilen Kommunikationsein- 
richtung mittels erster Kommunikationssignale, 
unter Verwendung eines nichtlinearen Kommunikationsmodells 
mit einem deterministischen Anteil und einem stochastischen 
10 Anteil, 

- welcher deterministische Anteil eine Abhangigkeit zwi- 
schen Kommunikationssignale einer Basisstation und Posi-^ 
tionen der mobilen Kommunikationseinrichtung beschreibt" 
und 

15 - welcher stochastische Anteil eine Unsicherheit des de- 

terministischen Anteils beschreibt, 
- bei dem unter Verwendung des zu der ausgewahlten Position 
zugehorigen ersten Kommunikationssignals der ersten Basis- 
station und dem nicht-linearen Kommunikationsmodell die 
20 ausgewahlte Position ermittelt wird, 



2.Verfahren nach Anspruch 1, 

bei dem mindestens eine zweite Basisstation, eingerichtet zu 
einer zweiten Kommunikation mit der mobilen Kommunikations- 
einrichtung mittels zweiter Kommunikationssignale, verwendet 
30 wird, wobei unter Verwendung des zu der ausgewahlten Position 
zugehorigen ersten Kommunikationssignals und/oder des zu der 
ausgewahlten Position zugehorigen zweiten Kommunikationssig- 
nals und dem nicht-linearen Kommunikationsmodells die ausge- 
wahlte Position ermittelt wird. 



wobei bei der Ermittlung eine nichtlineare Bayesianische 
Filtertechnik angewandt wird, wodurch ein mogliches Auf- 
enthaltsgebiet der ausgewahlten Position bestimmt wird, 
aus welchem die ausgewahlte Position ermittelt wird. 
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3. Verfahren nach An'spruch 1 oder 2, 

bei dem unter Verwendung des ersten und/oder des zweiten Kom- 
munikationssignals ein und/oder mehrere entf ernungsabhangige 
Parameter ermittelt werden bzw. wird, welche bzw. welcher von 
einer Entfernung der mobilen Kommunikationseinrichtung zu ei- 
ner von den Basisstationen abhangen bzw. abhangt und unter 
Verwendung des sen oder derer das mdgliche Auf enthaltsgebiet 
ermittelt werden bzw. wird. 

4. Verfahren nach einem der vorangehenden Anspruche, 

bei dem die Unsicherheit des deterministischen Anteils eine 
Unsicherheit der Kommunikationssignale, insbesondere ein 
Messrauschen, und/oder eine Unsicherheit der Abhangigkeit, 
insbesondere eine Unsicherheit des Kommunikationsmodells, 
sind bzw. ist. 

5. Verfahren nach einem der vorangehenden Anspruche, 

bei dem das zu der ausgewahlten Position zugehSrige erste 
Kommunikationssignal zeitabhangig ist und/oder fur einen 
Zeitpunkt k gemessen wird. 

6. Verfahren nach einem der vorangehenden Anspruche, 

bei dem die nichtlineare Bayesianische Filtertechnik iterativ 
angewandt wird. 

7. Verfahren nach einem der vorangehenden Anspruche/ 

bei dem als die nichtlineare Bayesianische Filtertechnik ein 
Gausscher Mischf ilter-Algorithmus oder ein Kalman Filter, 
insbesondere ein Extended Kalman Filter, verwendet wird. 

S.Verfahren nach Anspruch 1, 

bei dem der Gaussche Mischf ilter-Algorithmus ein Prior Densi- 
ty Splitting Mixture Estimator (PDSME) ist. 

9.Verfahren nach einem der vorangehenden Anspruche, 
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bei dem ein Benutzermodell verwendet wird, welches eine Bewe- 
gung der mobilen Kommunikationseinrichtung beschreibt. 

10. Verfahren nach Anspruch 9, 

5 bei dem unter Verwendung des Benutzermodells die Bewegung der 
mobilen Kommunikationseinrichtung fur einen Zeitschritt be- 
grenzt wird. 

11. Verfahren nach Anspruch 9 oder 10, 

10 bei dem unter Verwendung des Benutzermodells eine Posit ions- 
schatzung fttr die ausgewahlte Position ermittelt wird. 



12 .Verfahren nach Anspruch 11, 




bei dem die Positionsschatzung unter Anwendung der nichtline- 
15 aren Bayesianischen Filtertechnik ermittelt wird. 

13. Verfahren nach Anspruch 11 oder 12, 

bei dem die Positionsschatzung als ein Startwert fttr die Er- 
mittlung des mogliches Auf enthaltsgebiet der ausgewahlten Po- 
20 sition durch die Anwendung der nichtlinearen Bayesianischen 
Filtertechnik verwendet wird. 

14. Verfahren nach einem der vorangehenden Ansprtiche, 



bei dem das Kommunikationsnetz ein WLAN-, GSM- oder DECT-Net 
25 ist. 




15. Verfahren nach einem der vorangehenden Ansprtiche, 

bei dem als ausgew&hlte Position der mobilen Kommunikations- 
einrichtung ein Kennwert des moglichen Auf enthaltsgebiet s, 
30 insbesondere ein Schwerpunkt oder ein Erwartungswert des mog- 
lichen Aufenthaltsgebiets, verwendet wird. 

16. Verfahren nach einem der vorangehenden Ansprtiche, 
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eingesetzt zur Lokalisierung eines Mobiltelef ons in einem di- 
git alen, zellularen Mobilfunknetz, insbesondere in einem GSM- 
Netz, 

wobei die mobile Kommunikationseinrichtung das Mobiltelef on, 
die erste Basisstation eine gesprachsftihrende Basisstation in 
dem Mobilfunknetz ist. 

17.Anordnung zur Ermittlung einer ausgewahlten Position einer 
mobilen Kommunikationseinrichtung in einem Kommunikationsnetz 
mit mindestens einer ersten Basisstation, eingerichtet zu ei- 
ner ersten Kommunikation mit der mobilen Kommunikationsein- 
richtung mittels erster Kommunikationssignale, 
unter Verwendung eines nichtlinearen Kommunikationsmodells 
mit einem deterministischen Anteil und einem stochastischen 
Anteil, 

- welcher deterministische Anteil eine Abhangigkeit zwi- 
schen Kommunikationssignale einer Basisstation und .Posi- 
tionen der mobilen Kommunikationseinrichtung beschreibt 
und 

- welcher stochastische Anteil eine Unsicherheit des de- 
terministischen Anteils beschreibt, 

- mit einer Positionsermittlungseinheit, die derart einge- 
richtet ist, dass unter Verwendung des zu der ausgewahlten 
Position zugehorigen ersten Kommunikationssignals der ers- 
ten Basisstation und>dem nicht-linearen Kommunikationsmo- 
dell die ausgewahlte Position ermittelbar ist, 

- wobei bei der Ermittlung eine nichtlineare Bayesianische 
Filtertechnik angewandt wird, wodurch ein mogliches Auf- 
enthaltsgebiet der ausgewahlten Position bestimmt wird, 
aus welchem die ausgewahlte Position ermittelt wird. 

18.Computerprogramm mit Programmcode-Mitteln, urn alle Schrit- 
te gemafi Anspruch 1 durchzuftihren, wenn das Programm auf ei- 
nem Computer ausgeftihrt wird. 
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19. Computerprogramm mit Prograinmcode-Mitteln gemafi Anspruch 
18, die auf einem computerlesbaren Datentrager gespeichert 
sind. 

20. Computerprogramm-Produkt mit auf einem maschinenlesbaren 
Trager gespeicherten Programmcode-Mitteln, urn alle Schritte 
gemSB Anspruch 1 durchzufUhren, wenn das Programm auf einem 
Computer ausgeftihrt wird. 
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Zusaramenf as sung 

Verfahren und Anordnung sowie Coraputerprogramm mit Programm- 
code-Mitteln und Computerprogramm-Produkt zur Ermittlung ei- 
ner ausgewahlten Position einer mobilen Kommunikationsein- 
richtung in einem Kommunikationsnetz 

Die Erfindung betrifft eine Ermittlung einer Position einer 
mobilen Kommunikationseinrichtung in einem Kommunikationsnetz 
(Lokalisierung) . 

Dazu werden aus Kommunikationssignalen einer Basisstation, 
welche mit der mobilen Kommunikationseinrichtung in dem Kom- 
munikationsnetz in einer Kommunikation steht, unter Verwen- 
dung eines nichtlinearen Kommunikationsmodell mittels einer 
nichtlinearen Bayesianischen Filterung ein mogliches Aufent- 
haltsgebiet ftir die mobile Kommunikationseinrichtung ermit- 
telt. Aus dieser wird dann die ausgewahlte Position der 'mobi- 
len Kommunikationseinrichtung bestimmt. 

• i 

■ i 
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